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Observational evidence is presented for nonlinear inter- 
action between mesoscale internal Kelvin waves at the tidal 
- u!t or the inertial - uji frequency and oscillations of syn- 
optic - Q frequency of the background coastal current of 
Japan/East Sea. Enhanced coastal currents at the sum - 
u!+ and dif - u- frequencies: u}± — ujt,i ± fi have proper- 
ties of propagating Kelvin waves suggesting permanent en- 
ergy exchange from the synoptic band to the mesoscale cij± 
band. The interaction may be responsible for the greater 
than predicted intensification, steepen and break of bound- 
ary trapped and equatorially trapped Kelvin waves, which 
can affect El Nino. The problem on the parametric inter- 
action of the nonlinear Kelvin wave at the frequency uj and 
the low-frequency narrow-band nose with representative fre- 
quency f2 ^ tj is investigated with the theory of nonlinear 
week dispersion waves. 



1. Introduction 

Internal Kelvin waves play a significant role in the dy- 
namics of the oceans. There are two basic types of the 
waves: equatorially trapped and boundary trapped. Kelvin 
waves propagating on the equatorial thermocline participate 
in the adjustment of the tropical ocean to wind stress forc- 
ing [P/w/ander, 1990]. Extensive data on Kelvin waves have 
been obtained recently, motivated in part by possible con- 
nection between the initial stages of El Nino and equatorial 
nonlinear Kelvin waves which may precede this event. Any 
relaxation or reversal of the steady trade winds (the east- 
erlies) results in the excitation of a Kelvin wave, which can 
affect El Nino. In other words Kelvin waves play a critical 
role in generating and sustaining of the Southern Oscillation 
[e.'g.,Fedorov, 2000]. 

In the dynamics of the coastal oceans Kelvin waves are 
used to interpret such phenomena as the instability of along- 
shore currents, the generation and variability of wind cur- 
rents on the shelf, and upwelling [Brink, 1991]. Fedorov 
and Melville [1995, 1996] considered Kelvin waves trapped 
boundaries and showed that waves could manifest nonlinear 
properties that is steepen and overturn or break. A wave of 
depression deepens the thermocline and breaks on the for- 
ward face of the wave. A broken wave of depression can form 
a jump [also called shocks or fronts: Lighthill, 1978; Philan- 
der, 1990]. In turn, the breaking of Kelvin waves may be im- 
portant in mixing, momentum and energy transfer in coastal 
oceans. Other example of nonlinear Kelvin wave dynamics 
includes the problem of nonlinear geostrophic adjustment in 
the presence of a boundary [Helfnch et al., 1999; Reznik and 
Gnmshaw, 2002]. 
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The spectrum represents one of the major characteristics 
of the waves. It is used as a representative statistical descrip- 
tion of the wave field in studies of nonlinear interaction [e.g., 
Hibiya et al, 1998], acoustic propagation [e.g., Colosi et al., 
1998], and mixing parametrization [Polzm, 1995]. Filonov 
and Novotryasov [2005, 2007] studded the wave band of 
temperature fiuctuation spectra in the coastal zone of Pa- 
cific Ocean and observed that in the high-frequency band of 
temperature spectra the spectral exponent tends to ~ u)~^ 
at the time of spring tide and on the western shelf of the 
Japan/East Sea, in the uJi <^ uj ^ N, range, where A^* is 
the representative buoyancy frequency and cui is the iner- 
tial frequency, the spectral exponent tends to ~ uj~''. These 
features Filonov and Novotryasov [2007] simulated by the 
model spectrum of nonlinear internal waves in the shallow 
water. They considered interaction of high-frequency waves 
with the wave at the tidal frequency and shown that the 
spectrum of high-frequency internal waves take the univer- 
sal form and the spectral exponent tends to ~ uj~^. 

In this paper, we present observational evidence for non- 
linear interaction between synoptic oscillations of the back- 
ground current at the representative frequency Q and non- 
linear Kelvin waves at the tidal - u>t or the inertial - uJi 
frequencies. Findings are based on a well defined spec- 
tral peaks at the sum - u}+ and dif - oj- frequencies of 
inertial, semidiurnal and synoptic motions measured by 
current meters and temperature records collected in the 
coastal Japan/East Sea. We made rotary spectral analyses 
of these records and found that the clockwise rotary spec- 
trum for coastal currents measured at 35-m depth in the 
1999 year has well-defined spectral peaks at the frequen- 
cies ujt ~ l/12.4(cph), r2i,2 ~ 1/64, l/102(cph), as weU as 
sym and difference frequencies ijj± — ujs i: r2i,2 and their 
overtones uj„ — nuj±{n = 1,2,3). Analyses of coastal tem- 
perature records collected in the 2004 year showed that the 
spectrum of temperature variations has an analog form with 
spectral peaks at the inertial frequency uji ~ 1/17.8 (cph) 
and ill, 2 ~ 1/80, 1/160 (cph) synoptic frequencies, as well as 
sym and difference frequencies uj± — ctJ^ ±fii.2 and overtones 
LOn = nLj±{n — 1, 2, 3). 

With the theory of nonlinear interaction among week dis- 
persion waves [Gurbatov et al., 1990] we consider the prob- 
lem on the parametric interaction of the nonlinear Kelvin 
wave at the frequency uu and the low-frequency narrow-band 
nose with representative frequency i} <^ cj. We show that 
nonlinear interaction between Kelvin wave and no stationary 
low-frequency component of the background coastal current 
excited by the atmospheric forcing leading to the intensifi- 
cation of tidal and inertial currents on the sub-surfers and 
near-bottom layers of the coastal zone. Analogy phenomena 
of the Kelvin wave intensification can occur in the equatorial 
ocean. 



2. Observations 

In the first experiment, observations of internal waves 
were performed at the coast of the Gamov Peninsula area 
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Figure 1. (a) Study area on the Japan/East Sea shelf, 
September 1999, 2004. The mooring location in 2004 is 
shown in Arabic. The circle indicates the mooring loca- 
tion in 1999 and the location of the mooring vessel, from 
which hourly casts were conducted on September, 1999, 
2004. (b) Daily mean vertical profiles of temperature T 
and buoyancy frequency N . (c) Records of the meridional 
(light curve) and zonal (heavy curve) components of the 
current at the 35-m depth in autumn 1999. (d) Record 
of the temperature variations at the mooring buoy 1 in 
autumn 2004. 



The obtained data were analyzed by unified standard 
spectral techniques [Emery and Thomson, 1997]. The ap- 
proach involved (i) the elimination of low-frequency compo- 
nents with a Tukey's cosine filter, (ii) the splitting of the 
resulting series into nine 37.2-hour segments (three semidi- 
urnal tidal periods each), (iii) the calculation and averaging 
of spectral densities by segments, and (iv) the smoothing of 
the averaged spectral densities by a five-point Tukey's fil- 
ter. The number of degrees of freedom in the processing 
amounted to approximately 10 for the first experiment and 
20 for the second, providing a reasonable reliability of the 
results of spectral analysis. 

3. The model 

A spectrum model of nonlinear interactions among inter- 
nal Kelvin waves is developed. It is assumed that H/\ <^ 1, 
and A/H <^ 1, where H is the water depth, A is a character- 
istic wave-length, A is a representative wave amplitude. The 
basic component of this model is the simple wave equation. 
For the first most powerful mode of the small-amplitude 
Kelvin waves the equation is written as 



du/dx — audu/dr — 0, 



(1) 



where u{t, x) is the alongshore current, a:: is a horizontal co- 
ordinate, T = t — x/c, t is time. The parameters a and 
c are the coefficient of nonlinearity and the phase speed of 
long internal waves, respectively. Parameter of nonlinearity 
is determined by the background density and is related in 
the Boussinesq approximation as: 



of the Japan/East Sea in 1999 year. The records of current 
meter collected from a mooring deployed water depth 40 m, 
during 2 weeks in September (Figure la). To gauge temper- 
ature, speed and direction of currents, the buoy was instru- 
mented with the Russian-made POTOK integral instrument 
at a depth of 35 meters. The POTOK had a temperature 
measurement resolution of 0.05°C. The sampling rate was 
15 minutes. The time series of the meridional (light curve) 
and zonal (heavy curve) current components from this moor- 
ing are shown on Figure Ic. The Canadian Guideline CTD 
profiler was used for the vertical profiling of water. The 
hydrostatic pressure, temperature and electric conductivity 
of seawater were measured during profiling. The profiling 
errors were no greater than 0.01°C in temperature and 0.02 
psu in salinity. Figure lb shows the mean temperature and 
buoyancy frequency profiles for the study site. Details of 
this experiment (results, methods of their processing, etc.) 
were reported by Novotraysov, et. al [2003]. 

In the second experiment, measurements of internal wave 
band temperature fluctuations were performed during 18 
days starting 3 September 2004. The time records of tem- 
perature were collected from two moorings deployed along 
the coastline at a distance of 800m from it, approximately at 
40m depth and separated by a distance of 5.5 km from each 
other. The first of them was deployed at 28m depth and the 
second one, at 35m depth (below the surface). The moor- 
ings were equipped with digital thermographs made by the 
Russian manufacturer. The devices had a measuring pre- 
cision of 0.05° C for temperature. The sampling rate was 1 
minutes. The temperature and salinity vertical profiles were 
performed on 20-21 September, from a vessel anchored be- 
tween the moorings with the Canadian Guideline CTD pro- 
filer. In total, 25 hourly casts were made. 
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where 2; is a vertical coordinate, positive upward. The phase 
speed c and the amplitude function of the wave mode 4>{z) 
are determined from the solution of the eigenvalue problem 



d^(j)/dz^ + c''^N'^(z) (t> = 0, 
with boundary condition 

0(-H) = 0(0) = 0, 
and with the normalization 



(3) 



(4) 



0max — 1. (5) 

Equation (1) is a basic model for study of Kelvin waves 
interaction. Let us consider interaction nonlinear Kelvin 
wave with frequency lvq (tidal cjt or inertial uii) and narrow- 
band synoptic noise (later SN) with frequency fl <^ uio using 
the spectrum model of nonlinear internal waves described in 
terms of the equation (1). Let the alongshore current u{x, i) 
at the boundary of the coastal area x = be the superposi- 
tion of internal wave with frequency ujq, amplitude ^0 and 
the noise ■&{t) with typical frequency 7 << ujq 



u{t,x = 0) = uo{t) = Aocos {ojot + ip) +-&{t), (6) 
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Figure 2. (a) Normalized clockwise rotary spectra of currents versus for current measured at 35-m depth in autumn 
2004 near the Gamov peninsula. The numerical digits over peaks are their periods (h). (b) The increased fragment of the 
spectrum in the surrounding of the semidiurnal frequency. Dark numerals 12.2; 6.2; 4.3 are carrier periods (h). 



where ip - is a random phase with uniform distribution in 
the interval [— 7r,+7r]. 

We confine our analysis to the wave evolution stage, 
which is characterized by condition x < x^,, where x* = 
{aAou}())^^ . In this stage the front Kelvin wave appears and 
it is not accompanied by generation of internal solitons. We 
introduce parameter d — x/x^, which determine the similar- 
ity between a Kelvin shock-wave and the wave and consider 
the case d < 1. For this case the spectral density of the wave 
field u is 



a{(jj;x) = I {cxp[— ia;axi?(f)] — l}e*"'dt - 



E 



2iTiujax 



i"Jn{yjd/Luo) 
2niujax 



{exp[inip — iujax'3{t)] — 1} x 

X oxp[i(ct; — nuJo)t]dt (7) 



From equation (7) we can afhrm that the first term is the 
spectral density of the synoptic oscillations and the second 
term is the spectral density of the Kelvin waves distorted by 
the synoptic oscillations. 

We consider the spectral composition of the Kelvin wave 
near of the harmonics with the number n of the carrier fre- 
quency Wo- Given that 7 ^ ojo we make replacement cu at 
ULUo in the exponential rate of the equation (7). Then we 
obtain the equation for the harmonic with number n 



Un {t, X) — An{x) COs[nLJ()t + n{ip{t) — UJ()f3xii)(t))] (8) 

and her amplitude depending on the parameter d = Aoaujox 
is 
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Figure 3. Normalized spectrum versus for record of the 
temperature variations at the mooring buoy 1 in autumn 
2004 near the Gamov peninsula. The numerical digits 
over peaks are their periods (h). 



and stochastic components: i9(t) — acos(fit) +&ns{'t)- Let 
the amplitude of the phase modulation related with regular 
component M = anClx is small, that is M <^ 1. Let the 
Kelvin waves have narrow-band spectrum, and the synoptic 
noise width 'y <^ loq, besides the noise, amplitude and phase 
of Kelvin waves have Gaussian distribution. In this case the 
correlation function of the harmonic with number n taking 
into account the equation (8) is 



An{x) — 2Jn{nd)Ao/nd. 



(9) 



From the equation (9), we can affirm that the amplitude 
of the harmonic is decreased with increasing of harmonic 
number ~ 1/n. 

Perform the more detailed spectral analysis of the nonlin- 
ear Kelvin waves, when the synoptic noise consists of regular 



Bn(T,x) = { — {cosncuoT -\ |cos(nc<Jo + Ot) + 



-|-cos(najo — n)r]}}cxp[— n D^r/2], 



(10) 



where D^{t) = D^{t) + {aoJoxY {a'^s — Bnsir)). Here 
D^ is the structure function of the Kelvin wave phase. 
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Table 1. The periods of the peaks of the temperature variations spectrum (digit numerals) and the calculated periods (light 
numerals) 
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Bns{T) =< "dnsit + T)-&ris{t) > is the Correlation function 
of the synoptic noise. We perform the Fourier transform of 
the correlation function _B„ (r, x) and obtain the formula for 
the spectrum of the harmonic 



. Al ~ {AnaanuJoxf 
S„ (ijj; X) = —S„ [cj - ncjo) H x 

x{S„ (cj - {nujo + n)) + Snico - {riLuo - ^))}, (H) 



her 1,2,3 and 4 harmonics. Our estimations show that the 
magnitudes of these spectral peaks interrelated with number 
harmonics by the approximate relation 



where 



Snii^;x) = — 



exp[— n D^(r)/2] cos(cijr)dr. (12) 



From equation (10), we can afflrm that spectrum of the 
wave have harmonics with frequencies ncuo as well as har- 
monics with sum and dif frequencies uj± — nujo ± ^ and am- 
plitudes increasing with growth of distance traversed by the 
wave. This means that at the coastal zone between Kelvin 
waves and synoptic noise takes place nonlinear interaction, 
which is accompanied by intensification Kelvin waves with 
side frequencies uj± — nuJo ± ^■ 

4. Discussion and Conclusion 

Let us turn to the records of the current velocity. Fig- 
ure Ic presents records smoothed by Tukey window with 
the width 1 h of the meridional (light curve) and zonal 
(heavy curve) components of the current in the coastal zone 
of Japan/East Sea in autumn 1999 year. It is of interest 
to note a low-frequency component of the variation in the 
record of the current velocity. Let us analyze circular rotary 
current component with the rotary spectral analysis [Emery 
and Thomson, 1997] . The spectrum of the clockwise rotary 
(CWR) component of the velocity is 

S—{UJ) ^ 0, 125 - [Suu + Ovv ~ ^Quv), 

where Su-a and Sw are the one-sided autospectra of the u 
and V Cartesian components of the velocity and Q„„ is the 
quadrature spectrum between the two components. These 
spectra were determined by unified statistical spectral tech- 
niques using the algorithms presented in [Emmery and Thom- 
son, 1997]. The spectra have N degrees of freedom, where 
N > 10. Figure 2a shows the CWR spectrum of currents 
versus log lj normalized by maximal value. It is of interest 
to note the groups of the well defined significant spectral 
peaks in the surrounding of the semidiurnal frequency and 
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(with error < 20%). 

There are two peaks with light and heavy tops at the rep- 
resentative frequencies f2i ~ 1/64 and ^2 ~ 1/102 (cph) in 
the low-frequency band of the spectra. 

Figure 2b presents the increased fragment of the spec- 
trum in the surrounding of the semidiurnal frequency. It is 
of interest to note a significant well defined spectral peak 
at the Ws = 1/12.2 (cph) frequency surrounded by side 
peaks with dark end light tops. Not hard to make sure 
that the frequencies of these picks are determined by for- 
mula u)± = ojs ± f2i,2 with error equals of experimental er- 
ror uncertainty. In particularly the strong spectral peak at 
the difference frequency ui- = 1/12.2 — 1/64 (cph) suggests 
that the motion in study region is dominated by not only 
the tide movements, but and the inertial movements. The 
pick at the frequency LOi = 1/17.7 (cph) is formed by the 
inertial movements. Thus the spectrum has fine structure 
in the band 1/14 - 1/16.5 (cph), which don't segregate with 
help our dates and methods of analyses. 

Let us consider the structure of the spectrum in surround- 
ings of the harmonics of semidiurnal frequency lOs (Fig- 
ure 2b). As in the past picks at the carrier frequencies: 
ajsi,ci;s2 are surrounded by side picks at the frequencies, 
whish are determined by formula u)± = u)s\.2 ± ^1,2. 

Thus spectral analysis of the clockwise rotary velocity 
shows that her spectrum is determined by oscillations with 
the tidal frequency uj^ ~ 1/12.2 (cph) and the frequencies 
f^i,2 ~ 1/72, 106 (cph) from synoptic band and also that 
spectrum has the fine structure in the neighborhood of the 
frequency ojs with side peaks at the frequencies, which are 
determined by formula ct;± = ojs ± ili,2 The similar struc- 
ture has spectrum in the neighborhood of the first and the 
second harmonics of the semidiurnal frequency 10 s ■ 

Then turn to the records of the temperature of the coastal 
water of Japan/East sea collected in autumn 2004 year. 
Figure 3 shows the averaged spectrum of the temperature 
records collected from two moorings and normalized by the 
maximum value of the spectrum. Spectra are calculated by 
unified standard spectral techniques [Emery and Thomson, 
1997]. The number of freedom degrees amounted to approx- 
imately 20. 

Of particular interest are the groups of the significant 
spectral peaks at the inertial frequency and her first har- 
monic and also the peaks in the surrounding of the semidi- 
urnal frequency lOs. 
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Arabic numerals over spectral peaks equal the periods 
(hour) of these peaks. These periods (h) are placed in the 
Table 1. Digit numerals equal the periods of the spectral 
peaks and light numerals equal the periods calculated by 
the formula l/cij± — l/uJa,i ± i/i^i,2- 

As indicated by Table 1 the periods of the spectral peaks 
and the periods calculated by the formula have approxi- 
mately equal values. Thus the spectral analysis of cur- 
rent and temperature meter records of the Japan/East Sea 
coastal water shows that the spectrum is determined by 
oscillations with frequencies of near inertial tOi r; 1/17.8 
(cph) and tidal u>s ~ 1/12.2 (cph) frequencies as well as 
the frequencies f2i_2 ~ 1/72, 106 (cph) from synoptic band. 
The analysis shows that the spectrum has the fine struc- 
ture in the neighborhood of the frequencies iJi^s with side 
picks at the frequencies, which are determined by formula 
ci;± = iJi^si^i,2- The similar structure has spectrum in the 
neighborhood of the first harmonics of the inertial frequency 

UJi. 

As a conclusion, we showed that among nonlinear internal 
Kelvin wave with frequency u> and the low-frequency narrow- 
band nose with representative frequency Q <^ lj occur the 
parametric interaction, i.e. among the low-frequency band 
and high-frequency band of internal Kelvin waves occur the 
energy exchange. Hence, as a result of nonlinear interaction 
the energy of internal Kelvin waves in the coastal and trop- 
ical oceans can grow at the expense of the low-frequency 
narrow-band nose energy. 
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